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Thermodynamics of Maximum Transition Entropy for Quantum Assemblies 

David M. Rogers, University of South Florida 

This work presents a general unifying theoretical framework for quantum non-equilibrium systems. 

It is based on a re-statement of the dynamical problem as one of inferring the distribution of collision 
events that move a system toward thermal equilibrium from an arbitrary starting distribution. Using 
a form based on maximum entropy for this transition distribution leads to a statistical description 
of open quantum systems with strong parallels to the conventional, maximum-entropy, equilibrium 
thermostatics. A precise form of the second law of thermodynamics can be stated for this dynamics 
at every time-point in a trajectory. Numerical results are presented for low-dimensional systems 
interacting with cavity fields. The dynamics and stationary state are compared to a reference 
model of a weakly coupled oscillator plus cavity supersystem thermostatted by periodic partial 
measurements. Despite the absence of an explicit cavity in the present model of open quantum 
dynamics, both the relaxation rates and stationary state properties closely match the reference. 
Additionally, the time-course of energy exchange and entropy increase is given throughout an entire 
measurement process for a single spin system. The results show the process to be capable of initially 
absorbing heat when starting from a superposition state, but not from an isotropic distribution. 

Based on these results, it is argued that logical inference in the presence of environmental noise is 
sufficient to resolve the paradox of wavefunction collapse. 


I. INTRODUCTION 


Few topics draw more interest and debate than the 
relation between quantum mechanics, the measurement 
problem, and information theory. The problem of sewing 
these together into a coherent description of a causal or 
probabilistic universe lies at the foundations of physical 
law. 

These long-standing topics are receiving renewed at¬ 
tention due to the increasing precision of experiments 
on trapped quantum systems [l|-|6j] and the corresponding 
ability to probe the nature of decoherence processes. [3- 
[§] Parallel developments in nonequilibrium statistical 
mechanics of open systems are building an increas¬ 
ingly powerful toolset for predicting their often non- 
intuitive behavior. [10 12] Important studied applications 
include spectroscopyjUj[l4j], quantum computing. fl5j re¬ 
action rate theory, jlq mixed quantum-classical molecu¬ 
lar dynamics fl7j. and superconducting nanocircuits, flgj 

Recently, a general unifying theoretical framework for 
non-equilibrium thermodynamics was proposed for classi¬ 
cal systems. [l9| Such a framework has been a long sought 
goal for the nonequilibrium community. j20i . I2I ] It takes 
the form of a maximum entropy distribution over each 
transition step between the current and future positions 
and momenta, 



x' ,p' ,t + dt 



= Z^6((x'-x)-^At) 

e ~6^Tj 2dt / 2a2 + dt:i: s^iJ^ dt / 2 D) 


The two constraints in the exponent are the average de¬ 
viation of the square of the action functional, and the 
average change in total energy. Their Lagrange multipli¬ 
ers are <t~ 2 /2, one fourth the inverse diffusion constant, 
and /3/2, half the inverse thermal energy. When the 
Hamiltonian is in the form, H = p 2 /2m + V(x) (j|^y = 

—dp/dt — Q^r), Eq. [T] delivers the usual Langevin pic¬ 


ture of molecular kinetics, where momentum updates are 

o 2 

given by dp = (F — ^-x)dt + adW . The kinetic partition 
function, Z^ in , is a moment generating function that con¬ 
tains all the information about force/flux relationships. 
It makes an appearance in both the nonlinear fluctuation- 
dissipation and entropy production relations. mill] it 
also provides a physical interpretation of the Lagrange 
multipliers (collision rate, temperature, etc.) as constant 
forces responsible for bringing about equilibrium through 
noisy energy exchange with external systems. Equilib¬ 
rium properties like temperature, pressure, or chemical 
potential are only defined for the central system through 
the action of an external thermostat able to exchange en¬ 
ergy, volume, or molecules - a consistent interpretation 
of temperature out of equilibrium. mm 

Importantly, this work established that dynamic pro¬ 
cesses can serve as the foundation for statistical mechan¬ 
ics. In contrast to the fluctuation theorems, it does 
not require Boltzmann-Gibbs equilibrium distributions 
(or even stationary states) for initial and final densities. 
Rather, the dynamics is determined directly from maxi¬ 
mum transition entropy. Apart from being the stationary 
state of the dynamics, the Boltzmann distribution plays 
no special role. The log-ratio of forward and backward 
(inference) probabilities, 
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determines the ‘amount of irreversibility’ for the entire 
process. It consists of changes in the information entropy 
of the system plus heat rejected to the environment. The 
reverse process is determined by attempting to infer the 
starting distribution from the final one - a process subtly 
different from time reversal. Most notably, since f3 is 
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a Lagrange multiplier favoring transitions with smaller 
energy exchange, A E, this changes sign during inference 
to find starting states with larger energy. 


Based on the success of this work, it is reasonable to 
ask whether these results can be extended to describe 
quantum dynamics. The result would be a theory that 
directly models the statistics of transitions between quan¬ 
tum states, deriving thermodynamics as a consequence. 
One major obstacle stands in the way of such an ex¬ 
tension. The quantum measurement problem tangles up 
position and momenta, obscuring the classical picture of 
phase space. This means that it is not possible to de¬ 
fine a single path or trajectory which a wavefunction un¬ 
dergoes to transition from one time point to the next. 
This problem manifests in a variety of ways, including 
very high frequency oscillations in the propagator for the 
Wigner distribution and wavefunction collapse in matrix 
mechanics. 


This paper presents a novel statement of a transition 
probability distribution over quantum states. It retains 
much of the structure of the classical nonequilibrium the¬ 
ory referred to above. The distribution reduces to ordi¬ 
nary quantum dynamics when a collision rate parameter 
is set to zero. At nonzero rate parameter, the transi¬ 
tion probability contains an exponential bias toward tra¬ 
jectories representing exchange of energy with an exter¬ 
nal thermostatic reservoir. Additionally, a normalization 
constant for each step of the dynamics acts as a moment 
generating function for the energy fluxes, and a measure 
of irreversibility over a single timestep can be directly 
interpreted in terms of heat and work. 


In what follows, we will prove the unique properties of 
this approach, and analyze several of its consequences. 
Section [II] provides logical consistency arguments that 
eliminate alternative possible transition probability as¬ 
signments. Section fll Cl tlien introduces a novel quantum 
Andersen thermostat acting on an explicit model of the 
environment in an optical cavity experiment. Section HTTl 
goes over non-standard mathematical notation and prop¬ 
erties of transition superoperators used here. Results in 
section m compare the simulation of our physically mo¬ 
tivated reference system (containing an explicit environ¬ 
ment) with the maximum transition entropy distribution 
(which contains no explicit environment). Section [IV B 31 
provides a detailed picture of energy and entropy balance 
during a wavefunction collapse process which could oc¬ 
cur in a noisy interferometry experiment. The conclusion 
discusses interpretation questions and future extensions 
and applications of the results. The appendix contains a 
detailed proof that the maximum transition entropy ap¬ 
proach reduces to the density-matrix master equation of 
Caldeira and Leggett in the higlr-temperature limit. 


II. QUANTUM TRANSITION PROCESSES 

Exact quantum dynamics follows unitary evolution, 

m = u(t) = 

or, for the density matrix p = J2Pa\4’a){' t Pa\', 
p(t) = 

The unitarity of I/ -1 = means that both angles and 
distances between alternate state vectors, ip a , remain un¬ 
altered during time evolution. This also leads the eigen¬ 
values of p to remain constant, so that entropy is a con¬ 
stant of the dynamics. 

A modification introduced by Jaynes[25[ was to con¬ 
sider multiple alternative, unitary time evolution opera¬ 
tions, 

p(t) = Y / PrUr(t)pmht)- (3) 

r 

The final state density matrix then becomes a mixture 
of outcomes due to the unknown Hamiltonian that was 
applied to evolve the system in time. Unfortunately, this 
mixing picture was left incomplete, since it could not 
account for the temperature of the external system, and 
thus to lead to the thermal distribution. 

For defining U r , we can represent the effect of the ex¬ 
ternal system as a collection of random impulsive forces. 
Each will cause the central system to undergo a transition 
(in the position representation) 

ip(x) —> e^ xAp ^ p(x). 

In the continuous limit, the index r goes over into A p. 
These transitions are the short-time limit for any cou¬ 
pling Hamiltonian linear in the coordinate, x. Of course, 
other forms of environmental noise (e.g. periodic or spa¬ 
tially localized functions) are possible. 

Allowing for a quantum source of random impulsive 
forces, the most general transition operator will be the 
superposition, 

pf = J /(r)/*(s)e-^+i™ pe £A-;| s fi dr ds 

Measuring the amount of momentum actually exchanged 
will collapse the combined system onto a definite value 
of r, so that we need only consider noisy transformations 
of the form (Eq. [3]), 

p' = J |/(r)| 2 e-^ (# -™ /t) pe^ ( ^"™ /t) dr. (4) 

Jaynes analyzed this propagator in detail for the case 
when P r (equivalently |/(r)| 2 ) does not depend on ip or 
p. j25| We will term this case the “unweighted” distribu¬ 
tion. One central result is the subjective H theorem, 

S(t) < S(t') < S(t) -J2 Pr logPr ’ ( 5 ) 
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For the increase in the von-Neumann entropy. jl2l [26j 
S(t) = — Tr [p(t) log/?(£)], of the density matrix, /?, prop¬ 
agated from time t to time t' using Eq. [3j The physical 
content of the theory is that the classical uncertainty in 
the random noise applied during time propagation serves 
as an upper bound to the increase in the quantum en¬ 
tropy of the system, S(t') — S(t) < — ^ r P r \°gP r . In 
addition, Jaynes showed that the unweighted process al¬ 
ways contains the uniform (infinite temperature) distri¬ 
bution as a stationary state. 

The first central question in the theory of open quan¬ 
tum systems is, “How do we assign probabilities to the 
collision events?” Fermi’s golden rule is the traditional 
answer, but depends on assuming a directionality and 
time-dependence for the external energy. In a real inter¬ 
action process, however, more collisions will occur when 
the central system is moving with larger velocity. In ad¬ 
dition, the energy of the external system will decrease as 
the collision proceeds. The system’s energy change due 
to the impulse will then depend on the system’s initial 
momentum and the external force correlation time. 

The classical nonequilibrium theory [2dj modeled these 
random collisions by a maximum entropy distribution for 
P(r) = |/(r)| 2 , conditional on the average impulse size, 

(zf ) / an< ^ the ener gy change, (E' — E) = (A pp) /m. 
Correlation (not considered in this work for clarity) was 
also added by constraining the autocorrelation function, 
Ap(t) Ap(t-r) \ 

At At /• 

This lead to the unnornralized distribution (Eq. |T] in 
an explicit form, where p = ( p(t + At) + p(t))/2), 

P(Ap) oc exp Ap 2 /(2a 2 At) — ^App/m|. 

This result is easily shown to be equivalent to the 
Langevin equation: 


scheme for moving from a starting pure wavefunction to 
a density matrix a short instant later. The first approach 
leads to transition weights, (u r ip\H\U r ip\ — (ip\H\ip^, 
that depend exponentially on ip. Sec. Ill Al shows this 
choice to be inadmissible based on the known properties 
of the measurement process. The second method fixes 
this defect by working in superoperator space. Here, it is 
possible to average over all propagators, U r (-)U$, treat¬ 
ing them as members of a linear set of dimension N 2 
fEq. 1151) . Weighting the result by performing a similarity 
transformation is equivalent to weighting each propaga¬ 
tor, and results in a propagator that is linear in p except 
for normalization. 

A. Maximum Entropy conditional on ip 

Since each possible time evolution operation ( U r ) influ¬ 
ences the energy of the isolated system, we can choose P r 
in the quantum evolution operation (Eq. [3]) to be consis¬ 
tent with energy conservation. For a prescribed average 
energy change, the maximum entropy form is, 

P r OC p0 e -|<V'(0)|!7 r t Wr-il|v-(0))_ ^ 

The probability distribution over r is therefore condi¬ 
tional on the starting state, ^(0). This conditioning re¬ 
duces the likelihood of interactions with the external en¬ 
vironment that raise the energy of the system. 

Using random impulsive U r with this probability 
assignment, the wavefunction ip continually diffuses 
through Hilbert space. Despite the similarity to quan¬ 
tum state diffusion, [23-(29| it is not clear whether any set 
of evolution operators and constraints will bring the two 
into the same form. 

The analysis of Eq. [6] and the associated “kinetic” par¬ 
tition function, 


p’ — p = Fix) At + A p = (Fix) — ^— p)At + adW. 

2m 

In the quantum case, the dynamics are realized in terms 
of density matrices rather than operators. To find an 
analogous result, the next subsections develop this anal¬ 
ysis for several possible exponential biasing forms. 

Once a probability scheme is chosen, one encounters 
the second central question of open quantum systems, 
“Should time evolution be described by statistics of mul¬ 
tiple individual wavefunctions or by the statistics of the 
time-evolution operator?” This work refrains from con¬ 
sidering statistics among alternate density matrices or 
among non-orthogonal sets of wavefunctions. 

In principle, there are two different types of prior infor¬ 
mation for assigning probabilities to alternate transitions 
- creating weights that rely on complete knowledge of the 
system wavefunction, ip, or incorporation of a weighting 
into the transition operators themselves. These are dual 
to one another in the same sense that the Schrodinger 
picture is dual to the Heisenberg one. Both lead to a 


z kin (\,<p) = Y.p? e ~ XMAHW ’ 

r 


are relatively straightforward, except that the probability 
distribution is over the random perturbation rather than 
over the physical system of interest. 

Since the propagator will preferentially weight certain 
outcomes over others, it implicitly gives a model for wave- 
function collapse. Measurement happens when one of the 
set of possible outcomes at time t + At is selected by the 
environmental noise process. 

To be consistent with observation, however, the col¬ 
lapse should be completely described by the density ma¬ 
trix at time t. Extra information is not allowed because 
individual wavefunctions do not represent mutually ex¬ 
clusive events. [30| Rather, any two wavefunctions, say ip 
with pi and <p with p 2 , span a 2-dimensional subspace 
where the probability of any measurement can be repre¬ 
sented by arbitrary linear combinations [25| 


[• \f¥i = [•v/pTV' I Vpifi] 


cos 9 
sin# 


— sin 9 
cos 9 
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However, the choice of Eq. [G] leads to a change in the 
final density matrix ( p') when this reparameterization is 
done. This dependence would make it experimentally 
possible to distinguish different combinations of wave- 
functions making up the density matrix (i.e. 9). How¬ 
ever, this is not observed in practice. Therefore ip is un¬ 
suitable for use as prior information to directly determine 
the distribution of U r . 

We note that quantum state diffusion does not seem 
to have this drawback, since the form of the dynamics is 
such that the density matrix found by averaging \ip)(ip\ at 
each timestep conforms to the master equation in Lind- 
blad form, [la, Hl| 

In addition to this problem, note that there is also a 
loss of compactness in describing a weighting over tran¬ 
sition processes, U r . For example, two distinct propaga¬ 
tors, U r and U s , could lead to the same final ip' = U r<s ip. 
It would then be most economical to collapse these two 
propagators into a single logical possibility, reducing the 
upper bound of the subjective H theorem ([5]). This can 
be done using some linear algebra to state the transition 
probability in a way that does not depend on the (arbi¬ 
trary) space of random propagators, r. The bare mini¬ 
mum description of the transition probability for a wave- 
function is given by the 4-index superoperator [ljl HH 

(Sec. EH). 


G = J2 p rU r (-)U}. ( 7 ) 

r 

Expressions built in terms of this operator will then de¬ 
pend only on the actual action of the propagators on the 
system, rather than the abstract space of r. 


B. Energy-Weighted Propagator 

This section examines a probability assignment based 
directly on a weighting of the superoperator, G, toward 
transitions representing larger amounts of energy loss 
from the system. First, re-express the transformation 
(Eq. Q in the energy basis as p\,-, = V,, ( ; iy . 
Next, apply an exponential bias toward propagators with 
larger amounts of net energy loss. This leads to a re¬ 
weighting, 


G,/ 




= G° , p ~H^E i +AE j )/2 
i'i \jj’ 


( 8 ) 


This is equivalent to transforming the N 2 canonical prop¬ 
agators derived by an eigenvalue decomposition of G°. 
For proofs and related notation, see Sec. 1IIII Finally, we 
arrive at a prescription for adding a bias accounting for 
energy conservation to an initial stochastic transition su¬ 
peroperator, G°, 

G = ( e - A ^(.) e Ai7/2) G 0( e AiI/2(.) e -Aif/2). (g) 

The description of the dynamics then becomes some¬ 
thing like 


P = G[p\. 


However, p' is not necessarily normalized. Instead, the 
overall normalization can be expressed as a matrix oper¬ 
ator, 

Tr [G[p\] = TVE G i'i\n'Pij\ = Tr \ Z P\ 
ij 

Z = (Tiq_[G]) T = J2P?e XH/2 U}.e- XH U r e XH/2 . 

r 

( 10 ) 

In principle, there are three options for normalizing 
p'. The simplest is applying a single scaling constant to 
the whole input density matrix, G$[p] = G[p]/ Tr [Zp]. 
Failing that, the density matrix could be scaled before 
applying G to give a linear superoperator, 

G ( $[p] = G[Z- 1 / 2 pZ- 1 / 2 }. (11) 

As a final alternative, the starting density matrix can be 
normalized by decomposition into a set of pure states, 
p = Y2 k Pk\'4’ k )('ip k \, each of which is given an individual 
normalization constant, 

g n = mz\i>k) 

k 

Gn[p] = '^2pkG ( ^[\'ip k ){'tpk\\- (12) 

k 

The first choice has nice mathematical properties, but 
is physically inadmissible. It shows non-causal behavior 
by combining otherwise independent parts of a statistical 
mixture of states. Specifically, the single normalization 
constant takes contributions from all starting states in a 
mixture, 


The weights connecting transitions between energy eigen¬ 
states are of maximum entropy form. It is difficult 
to prove a more general statement about the entropy 
over transition superoperators. Since a more general 
extremum principle has not yet been found, this ap¬ 
proach will be termed the energy-weighted, propagator in 
the analysis below. 

When written in an arbitrary basis, the bias on the en¬ 
ergy change (Eq. [ 8 J has the effect of transforming each of 
the unitary operators of Eq.[3]into U r —> e~ XH GjJ r e XH ^ 2 . 


This operator does not make a suitable transition dis¬ 
tribution, since summing over final states does not give 
back the starting probabilities, p k ■ Consequently, it can¬ 
not reproduce the starting density matrix, and will not 
be considered further. 

The second choice defines a linear superoperator that 
makes a more suitable transition distribution, but does 
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not act like a moment generating function as we would 
expect from maximum-entropy statistical mechanics. In 
particular, the first derivative of the matrix logarithm 
does not give the energy change over a time-step unless 
Z has a special form (commutes either with the energy, 
H , or the density matrix, p). 


Tr 


d log Z 
<9A P 


E p ° Tr 


[-ff,t/ r ]Z- 1/2 pZ _1/2 E/? 


This 

tive, 


equation is proved by noting the deriva- 


STr 


G\ [ p] /d\ is zero and using Z - 1 / 2 = 


exp(— \ logZ). Also, the starting density matrix is not 
recovered using the naive definition, 


[G^PpU w = G v 


7 -l/2 n z -l/2 

,’n\mj ,Zj ni i 


That leads instead to Tr<_[G)y p)ij = SijPij. Although 
it will not be considered in this work, numerical experi¬ 
ments (not shown) find that this operator behaves much 
like Gjv, when the noise level in P® is much smaller. 

The third choice propagates orthogonal wavefunctions 
belonging to the initial density matrix separately. This 
leads to a proper causal picture[25| where a single pure 
state chosen from a statistical mixture at each point in 
time is responsible for the system’s future behavior. It 


( 2 ) 

reduces to the second choice, G N , when Z commutes 
with p. 

It is also simple to check that propagating the mixture 
of Eq. [7] gives an answer independent from 9, 


p'= v \G N [W)m\+ P ' 2 G N [W)^\] 

= (Pi l(V#')| 2 +P 2 I<V’I0 , )I 2 ) G! jv[IV’)(V’|] 

+p'2\{mt)GNum\ 

= PiG N [\ip){i/j\] +P2G4 1 </>}(</>!]. 


Interestingly, this means that, although G must be spe¬ 
cialized for p , the choice of propagating only the eigen¬ 
states of p in Eq. [12] is not unique. The same final den¬ 
sity matrix, p', would result from propagating any non- 
orthogonal weighted set of starting states which combine 
to give p (an array in the terminology of Ref. El). 


1. Transition Free Energy and Entropy 

The remainder of the paper introduces a reference pro¬ 
cess and proceeds to analyze the qualitative and analyt¬ 
ical mathematical properties of the third choice, Gjv • In 
particular, we find that a kinetic free energy functional, 

F = -^Pklog^klZli/jk), (13) 

k 

(where Z was defined in Eq. 1101) has properties consistent 
with maximum entropy thermodynamics. 


A nonnegative relative information entropy can be de¬ 
fined for each transition, 

dStot = Tr [G N p (log G N p - log p'G' N )\. (14) 

By the Klein inequality, this is necessarily positive when 
p'G' N is any normalized, completely positive 4-index ten¬ 
sor. 

Conventionally, G' N is chosen to be a “reverse” trans¬ 
formation from the final density p' to an inferred joint 
density matrix - so that dS'tot is a measure of the ir¬ 
reversibility of the process during time At.JH, [IsJ This 
choice appeals to information theory on the grounds that, 
after discarding p in favor of the updated p ', G' N presents 
a “best guess” as to the 2-time probability density of the 
system. [191] Since dS’tot is the information lost during the 
update step, it must obviously go to zero for deterministic 
evolution and depend on the ratio of forward to reverse 
probabilities. 

Physically, Eq. [T|] should be derived by considering 
the trade-off between spontaneous decrease in the sys¬ 
tem’s information entropy and heat released into the 
environment. 24] In this case, the reason for using a re¬ 


verse transformation is so that G' depends on temper¬ 
ature like e^ A£; / 2 . The difference, logG — logG', then 
contains a heat release term (negative of the system’s 
energy gain (/3AE)) that counter-balances increase or de¬ 
crease of the system’s instantaneous information entropy, 
AS’inf— (/3AE) ~ AS'tot > 0. This should be compared to 
the classical result (Eq. [2]), where, for quantum systems, 


Sinf = - Tr [p log p\ 


Re-stated in terms of the system itself, the system cannot 
store thermal energy without increasing its information 
entropy. 

A straightforward choice based on these considerations 
is, 

G' = (e A/2p (.)e" A/2p )G°(e" A/2p (.)e A/2p ), 


and consequently 

G)v = E( I^X^I(-)/)G , (/(-)|V’l.)<V’l.|)/<^|G'[/]|^). 

fc 


Note that the direction of time plays no special role other 
than swapping the application order of p' - matching the 
final indices i',j' so it encounters U\-)U. The results 
in Sec. IIV B 31 show that this definition is most closely 
related to the information loss interpretation, since it de¬ 
composes approximately as, 


dStot — (IS] ,, f + f3dQ irrev- 


C. Reference Process 

To illustrate the features of the energy-weighted prop¬ 
agator, it is helpful to compare against a reference pro¬ 
cess. Many recent experiments on molecular beams and 
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trapped atomic systems have been carried out to probe 
tire dynamic properties of wavefunction collapse. [3?, 0, 
E, EH Although experiment makes the best reference in 
principle, direct comparison is limited by uncertainties 
in the source and strength of environmental noise. These 
differences, reflected in alternative choices for U r , make 
noticeable differences as the noise strength increases. 0 
In addition, the Hamiltonians for most experimental se¬ 
tups introduce additional complexities which are not es¬ 
sential for understanding the theory. 

We therefore consider a physically motivated theoret¬ 
ical reference process for environment-induced decoher¬ 
ence and dissipation. This provides a strong intuition for 
the role of energy-weighting in the stochastic propagator 
of the density matrix. 

The simplified reference system consists of two reso¬ 
nant harmonic oscillators. [35j Oscillator a represents the 
system of interest (e.g. the atom in a Paul trap). Os¬ 
cillator b represents a single electromagnetic mode of an 
enclosing cavity. Spin is ignored, and the two are simply 
coupled by an energy term, cx a Xb- The Hamiltonian is 
thus: 


a) 


c) 


a 

0 ) 

m 

m 

s 



02468 02468 


^cavity 


H = hjjjicJa + V*b + 1) + c--(a* + + b), 

2 raw 

where d, b are the lowering operators for the uncou¬ 
pled harmonic oscillators. Similar models using a har¬ 
monic oscillator as the environment and a 2-level atom 
as the system of interest have also been investigated in 
the context of deriving fluctuation theorems for exact 
dynamics. [H, E3| Further elaborations on this model to 
include the atomic spin and the direction and polariza¬ 
tion of the field |38 h 4C)I ] could be used for comparison to 
experiment. @ 

The environment in our reference model has the re¬ 
sponsibility for introducing both decoherence and energy 
dissipation. To model this, we propose a physically mo¬ 
tivated thermostatting scheme. It is based on the Ander¬ 
sen thermostat for molecular dynamics. There, the ve¬ 
locity of a particle is periodically replaced with a random 
sample from the Boltzmann distribution. (4l| By analogy, 
this scheme periodically replaces the cavity state with a 
pure energy level sampled from the Boltzmann distribu¬ 
tion. 

The process can be visualized as in Fig. [L] The 
joint wavefunction is represented in the uncoupled en¬ 
ergy space as 

\4>)(t) = ^2 ' t l’nm{t)\n)a®\ni)b- 

n,m—0 

The panels of Fig. [L] show probability densities for these 
joint states, \ipnm\ 2 - Starting from a state where both 
energy levels are known in[T^, the joint system undergoes 
deterministic evolution according to 

I ^)(t) =e~h H \^){0), 


FIG. 1. Schematic of the quantum Andersen thermostat (ref¬ 
erence process), (a-c) Exact evolution of the coupled sys- 
tem+cavity from a Foc.k state quickly leads to an entangled 
state, (c) Measuring the cavity energy selects a subsample 
of the system, removing coherences.(d) Replacing the cavity 
state with a thermal sample adds or removes energy. The 
thermal nature of the environment is responsible for dissipa¬ 
tion. 


which introduces a system-field superposition (P>c). 
Measurement of the cavity mode JTJ;) collapses the sys¬ 
tem into one of the pure cavity states, ip a <S> | m)b = 
C^2 n V , nm(t)jn) a ) ® | m)b, chosen in the usual way with 
probability P(m) = J2 n l^wl 2 ■ Finally, the cavity mode 
is replaced (jTJi) with a new state sampled from the ther¬ 
mal distribution, 


ipa ® \m)b -t V’a <S> \m')b, P{m!) oc e /3RaW . 


There is a distinction in this process between deco¬ 
herence and dissipation. Decoherence occurs when the 
environmental state is measured, since a state containing 
all the information about the the system-cavity super¬ 
position is projected onto a single cavity mode. This 
step usually introduces small increases in the energy of 
the joint system because the coupling, cx a Xb has only 
off-diagonal terms in m - which are removed by the pro¬ 
jection. However, replacement of the cavity mode with 
a thermal sample accomplishes energy exchange heat¬ 
ing or cooling the system until thermal equilibrium is 
obtained. 
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1. Deterministic Limit 


so that 


The measurement process described above incessantly 
projects the environment into a pure state. This type 
of measurement traces over the environmental degrees of 
freedom, so that at each measurement point the environ¬ 
ment is seen to follow a deterministic path without su¬ 
perpositions. Starting with random off-diagonal matrix 
elements for the environment will have the same effect - 
canceling paths that begin in a superposition state. 

If the trace (measurement) process is made continuous, 
the environment follows a definite trajectory, and is only 
felt as a deterministic external force - replacing H a b with 
H a + V(x a ,Xb)- This is the quasi-classical model used 
for many time-dependent quantum processes, where the 
Hamiltonian depends predictably on time through Xb- 


III. METHODS - SUPEROPERATOR 
NOTATION 


For analyzing the re-weighting of Eq. 0 we introduce 
the following notation for density matrix superoperators. 
The canonical form for a superoperator is[26| 

0(N 2 ) 

G= E 7 kC k (-)Cl (15) 

fc=l 


It is a 4-index quantity. Explicitly, G^i ny = 
Yhk'YkGk,i'iGkj'j- Unless otherwise noted, each index 
(i,j,etc.) ranges over integer energy levels, 0, • • • ,N — 1. 
Superoperators appropriate for the density matrix can all 
be written in the form of Eq. [15] by virtue of the (super- 
Hermitian) symmetry property Gj/qy/ = G*,^,. This 
property also guarantees a spectral decomposition into 
N 2 matrices, {Gfc}, which are mutually orthogonal - as 


defined by J2i>i Gk,i’iC* i , i = TV 


c}c k 


= 6 , 


kl • 


We define the inner and outer trace operations, 


T 

ij 

TV[G]^=£T iW . 

i’j’ 


Superoperators are applied to matrices by a substitu¬ 
tion pattern, 


(A(-)B)[p]=ApB, 

and are composed following the same pattern, 

(A(.)B)o(C(.)D) = (A(.)B)[(C(.)D)[p}] 

= (. AC{-)DB)\p\. 

We also define a simple algebraic notation for multipli¬ 
cation, 


(A(.)B)(C(-)D) = (AC(-)BD) 


(A(-)B) o (C(-)D) = ( A(-)I ) (C(-)D) (I(.)B ). 

Note the difference in application order between com¬ 
position and multiplication. The multiplication opera¬ 
tion makes certain manipulations clearer. In this nota¬ 
tion, the Choi-Jamiolkowski isomorphism is simply ob¬ 
tained by writing the identity for a complete basis, 

G='£mm\)G(\i)mi). 

ij 


Writing an initial transition tensor explicitly in the en- 
ergy basis, {|fc)}^ -1 , 

G °= E (l* / )(*l(‘)li)0’ , l) > 

We define the direct product of a superoperator and a 
density matrix, p, as the projection onto the eigenspace 
of the density matrix ({pk, ipk}^ 1 ), 

Gp = E^( J (-)I^I)(^II) G (I^)(^I(-K), 

k 

so that G\p\ = Tr_y[Gp]. This operation effectively spe¬ 
cializes a superoperator to be consistent with its operand, 
p. The choice to project or not does not affect the dis¬ 
cussion in Sec. UT1 leading to Gat. 

The numerical results presented here made use of tabu¬ 
lated G® w tensors calculated by numerically integrat¬ 
ing Eq. U for |/(r)| 2 | = exp(—r 2 /2cr 2 )/ V2i rcr 2 using a 
20-point Gauss-Hermite quadrature rule and an exact 
matrix exponential computed through eigenvalue decom¬ 
position. 


IV. RESULTS 

The first role of the environment is to provide a mech¬ 
anism for decoherence by introducing random noise into 
the time evolution. This section presents results from the 
reference process, which show the importance of the to¬ 
tal energy in determining intermediate states just before 
collapse. 


A. Dynamics of the Reference Process 

The reference process makes quantum jumps through 
periodically projecting the cavity into a Fock state. Al¬ 
though this always leaves the central system in a super¬ 
position state, Fig. [2] shows that the probabilities for 
each system state still appear to follow a jump pro¬ 
cess. The resulting trajectory is very reminiscent of 
fluorescence measurements of single-atom excited state 
probabilities, 





We find numerically that the Andersen-type thermo¬ 
stat scheme employed for the reference process leaves the 
system very nearly in the equilibrium state of the coupled 
system. However, when the coupling strength is large, 
projecting out the cavity from the combined equilibrium 
density matrix gives a system density that can deviate 
significantly from the Boltzmann distribution (Fig. [I]). 
Practically, this means that (for physically correct rea¬ 
sons), the thermostat does not work at strong coupling 
and only yields the exact Boltzmann distribution in the 
zero-coupling limit. We shall find similar behavior for 
the maximum-entropy jump process. 

We first determine the dependence of the decorrelation 
time and energy dissipation in the reference process as a 
function of the coupling strength, c, and the time between 
interactions, i rese t- Fig. [2] shows the average system en¬ 
ergy level vs. time over 1000 stochastic trajectories, each 
starting from the ground state. The index labels the 
system energy level, projected from the complete sys¬ 
tem plus cavity state. Each parameter set was run for 
a length of 1500xf reS et, using N = 10 energy levels for 
each system and a numerical matrix exponent between 
each reset. [l3| The time-course of po fits well to an expo¬ 
nential with relaxation rate that increases approximately 

, 3/2 

^ reset' 

We have verified numerically that the relaxation rate 
depends on the square of the coupling strength, c, as pre¬ 
dicted by perturbation theory. [10] Finally, a dependence 
on temperature will come from the dependence of the 
jump probability on energy level. Since the jump proba¬ 
bility between n, m +1 and n+ 1 , m states by perturbation 
theory is proportional to nm, we might expect quadratic 
dependence on (m) + 1/2. However, at high temperature 
the jump probability empirically shows a linear, rather 
than quadratic dependence on ( n). Fitting the observed 
relaxation rates, r, using this scaling resulted in reason¬ 
able agreement for the constants, 


r ~ 0.0510747- 


1 / 2 ) 


(ustreset) ^ 


0.000154756. 

(16) 


B. Properties of the Energy-Weighted Propagator 

This section begins by showing that the qualitative 
dynamics and approach to equilibrium of the energy 
weighted propagator is very similar to the reference pro¬ 
cess above. Next, simple examples of the kinetic partition 
function are given. The section concludes by showing the 
relation between work, heat and entropy production. 

Fig. [3] shows the time-dependence of a single-harmonic 
oscillator simulated with the maximum-entropy matrix 
propagator derived in Sec. IIIBI The variance of the 
random impulsive forces was chosen to give the relax¬ 
ation rate corresponding to the reference process (us¬ 
ing Eq. m- The variance of momentum jumps is re¬ 
lated to the diffusion constant of Langevin dynamics as 


a 2 = 2D At. [24j As should be expected, the energy relax¬ 
ation rate was found to be linear in the diffusion constant, 
at r = 1.423971? per step. 

The dynamics of the weighted propagator appears 
qualitatively similar to the reference process. Fig. [3 Jj 
shows jump behavior that indicates the propagator in¬ 
duces localization of the wavefunction in energy space. 
The transitions are slightly more random, as might be 
expected from a maximum uncertainty model for inter¬ 
action with the environment. The average probabilities 
in panel [3 Jd (dashed lines) fit the reference process (solid 
lines) noticeably better than a simple exponential - com¬ 
pare to the green line in Fig. [2] We stress that the only 
fitting parameter is the diffusion constant, which was cho¬ 
sen so that the relaxation times agree between the two 
processes. 


1. Stationary Density Matrix 


The propagator G^/ admits the Boltzmann-Gibbs 
distribution with /3 = A as a stationary state of 
the dynamics. This is shown since G[e~ XH ] = 
e~ XH / 2 G 0 [I]e - XH / 2 = ce~ XH , and G°[/] = I. The un¬ 
weighted propagator always preserves the infinite tem¬ 
perature distribution, since it is an average over unitary 
operations (Eq. [3j). 

The propagator Gn (Eq. [12]) would not normally be 
expected to exhibit a Boltzmann-Gibbs stationary dis¬ 
tribution. However, empirically the stationary state of 
Gn at small A is very close to the Boltzmann-Gibbs dis¬ 
tribution with j3 = 2A (Fig. [4]). At large A, the stationary 
state becomes more like the Boltzmann-Gibbs distribu¬ 
tion with /3 = A for G^. 

The stationary distributions in Fig. [4] were found by 
repeated eigenvalue analysis on the re-normalized propa¬ 
gator, Gn- Convergence was achieved when the eigen¬ 
vectors of p rotated with N — I (V’felV’fc) I < 10~ 4 . 
The stationary state shows two regimes. The high- 
temperature, low diffusion constant limit gives a density 
matrix approximately proportional to e~ 2XH . The low- 
temperature, high diffusion constant limit gives instead 
the density matrix e~ XH . This behavior is evident at the 
break in the slope of probability vs. A. 

This behavior can be justified when there are enough 
energy levels in the system such that Gf m _ n)ra|m(m _ n) = 

G °(m+n)m\m( m +n .) > with G \m+n)m\m(m+n) decaying with 

increasing energy difference, n. In this case, the moment 
generating functions, 

W (Aj—\^G° i e -HEn-E m ) 
yrm(A) — / y tjr wrra | mrl e 


are nearly independent of the starting state, n, for small 
A. If they are independent, the Boltzmann-Gibbs distri¬ 
bution obeys the detailed balance condition, 


e ~2\E n 


qO 

mn\nm 


e 


-A (E m -E n ) 


= e ~ 2XE ' 


QO g-A (En-Em) 

nm|rari 


W n ( A) 


W m ( A) 
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FIG. 2. Panel (a) shows occupation probabilities for the ground and first excited state of a central system interacting with a 
cavity under quantum Andersen thermostatting (foj = 1, /3 = 0.1, c = 0.01, t rese t = 10). The pi curve is shifted down by 0.2 for 
visual clarity. Panel (b) shows the time-course of the average occupation probability for ground (po, upper curves) and first 
excited states (pi, lower curves). The parameters are the same as panel (a), except /3 = 1. Solid lines are averages over 1000 
trajectories. Dashed lines show exponential relaxation with a rate constant fit from po(t). 



FIG. 3. Dynamics of the weighted energy exchange propagator, (a) Single trajectory (/3 = 0.1) integrated by choosing a 
pure state from the density matrix at each time point. The pi(t) curve has been shifted down by 0.2 for visual clarity, (b) 
Comparison of combined density matrix evolution (smooth curves) with analogous quantum Andersen thermostat result from 
Fig. [lb- The temperature, (5 = 1, translates to A = 1/2 (Sec. IIVB 111 . 


so the equilibrium probability of energy level m is n m oc 

e -2XE m _ 

As the coupling strength (A) is increased, the mo¬ 
ment generating functions tend to their extreme value, 
W m { A) —> Gp m | m0 e A£ ' m_AE °. Detailed balance then gives 

7 r m oc e -A£ ' m-1 ° sG ° m i m0 . The probability of jumping to 
the lowest energy level, Go m i m o, depends on the collision 
rate parameter. As the rate parameter decreases, the 
cross-over point where the effective temperature of the 
stationary state doubles (from /3 = 2A —> A) is delayed 


to ever lower temperatures. This behavior is seen clearly 
in the numerical results. 

The upper left panel of Fig. [I] shows that an analo¬ 
gous effect happens in the reference model. There, the 
thermal state with the coupled system-cavity Hamilto¬ 
nian depends on the temperature and coupling strength. 
When the coupling strength is small, tracing over the en¬ 
vironment gives very nearly the Boltzmann-Gibbs den¬ 
sity matrix for the system. As the coupling strength in¬ 
creases, the system’s projected density matrix plateaus 
to a constant, minimum, temperature at smaller /3. This 
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FIG. 4. Temperature-dependence of the stationary states of (a,c) reference Hamiltonian and (b,d) energy-weighted propagator. 
Each color corresponds to a different coupling parameter (in a,c) or diffusion constant (in b,d). For (a) and (b), the upper set 
of lines show the first excited state log ratio, logpi/po, while the lower set show the second excited state ratio, logp 2 /po- The 
Boltzmann-Gibbs distribution would show straight lines with slope -1 and -2, respectively. Numerical noise at low temperature 
in (b) is due to difficulty converging the very low probability of the second excited state, p 2 - Panels (c) and (d) show the total 
magnitude of the off-diagonal elements of the 10 x 10 density matrix. 


may explain the requirement for a rescaled effective tem¬ 
perature in interpreting electron transport experiments 
probing the fluctuation theorems. fUf 


The lower panels of Fig. 0] show the total magnitude 
of all off-diagonal elements in the density matrix as a 
function of temperature. Although these are on the or¬ 
der of 1CT 3 (N = 10), both the reference process and 
the energy-weighted propagator show a maximum at in¬ 
termediate temperatures. In the reference system, the 
off-diagonal elements are mainly two above the diagonal, 
and can be traced to the double-excitation process, a'b 
in the linear coupling, x a Xb- 


2. Moments of the Flux Distribution 


The function F(— /3, p) (Eq.[13]) acts as a cumulant gen¬ 
erating function for the distribution of the energy flux. 
This is easy to show using successive differentiation of 
Tr [Zp\. From the definition of G (Eq. 0, 


d n Tr [Zp] 
d(-X) n 


Y — 


2 n i\j\k\l\ 


Tr 




G((-H) k (.)H l )[p\ 


E (") Tr [WGK-Hr^p]] . 
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FIG. 5. Energy change per unit time for a system initially 
in an energy eigenstate. The ground state (n = 0) absorbs a 
small amount of energy, and is plotted separately on the left 
scale. At small /?, other states (plotted on the right scale) 
initially absorb small amounts of energy but eventually lose 
energy quickly at lower temperature. Lines give the result 
from differentiating F (Eq. 11311 , while the points show a sep¬ 
arate calculation as the average energy change from a single 
time-step of dynamics. 


The summation set, £, contains nonnegative integers, 
i,j , k, l that sum to n. 

The cumulants of the energy flux (derivatives of F ) 
therefore break into sums over cumulants from each start¬ 
ing state, a, 


f) n F' 

= Em-i)"*-? (u) 

a 

While each cumulant, K™ is defined in the usual way 
in terms of n th order polynomial expressions in the mo¬ 
ments, 


K}= Tr HG n [ |VbXVb|] -(VbWi 
K? = Tr H 2 G N [\^j)(^\\ + Uj\H*\t 


- 2Tr 


HG n [H IVbXVb'l] ~(K)) 2 


etc. 


These are a suitable definitions for the n th order cumu¬ 
lants of the distribution of AH. na When p and p’ are 
both energy eigenstates, it reduces to the classical statis¬ 
tical distribution (45j . 


3. Wavefunction Collapse 

To illustrate the properties of the entropy, we give de¬ 
tails on the process of wavefunction collapse with the 
energy-weighted propagator. This will be illustrated with 
a simple molecular beam experiment (Fig. [6|).[46j Here a 
stream of atoms in a maximally entangled spin state are 
subjected to a random perturbing magnetic field during 


flight from a source to a detector. The final density ma¬ 
trix can be measured by counting the number of spin-up 
particles observed along each direction. 

A uniform external magnetic field provides the energy 
function needed to induce wavefunction collapse. Con¬ 
ventionally, the field is said to ‘measure’ the particle’s 
spin state. Because of the field, the particles arriving in 
the detector are only found in one of two states - up or 
down along the 2 -axis. This investigation of thermaliza- 
tion therefore demonstrates the hallmark of wavefunction 
collapse - that measurement forces the system onto the 
diagonal in the measurement basis. 

The mathematical properties of this system are simi¬ 
lar to recent experimental investigations on decoherence 
rate of quantum systems achieved using Ramsey interfer¬ 
ometry on trapped atoms in superpositions of coherent 
harmonic oscillator states. Q There, each of two states is 
tagged with a spin. A random displacement shifts the 
relative phase of the two states, resulting in decoherence 
that can be probed by moving the states back together 
spatially and then measuring the probability of up-spin 
along each detector angle, 9. 

The propagator corresponding to evolution of a spin in 
an external magnetic field is U r = e lr ' G , where a repre¬ 
sents a vector of Pauli spin matrices, and the direction 
is proportional to the magnetic field, r = B. The 

constants g s Pv stand for the appropriate gyromagnetic 
ratio and magneton. |47| Its operation can be computed 
exactly using a general formula for the exponential of a 
trace-free 2x2 matrix. (43] 

e lR = cos \R\I + i ^j 1 #, (18) 

\R\ 

where R = r ■ a = r x a x + r y a y + r z a z and \R\ = 
\j r x + r y + r z is the magnitude of both (positive and 

negative) eigenvalues of R. The action of the unitary op¬ 
erator can be calculated on a 2 x 2 Hermitian matrix by 
expanding as p = al + C, where C = c - c r, 

e iR (aI + C)e~ iR = p+ [e iR ,C]e~ iR (19) 

= p+ i si nWc°s|fl | [R , 01 

+ 2 !! j^pr (c-rfl-liiPC) ( 20 ) 

Its average effect on the density matrix is determined 
by the statistics of the direction. Integration over r 
is difficult in general, but we can derive useful results 
by choosing r x = r xy cos <j>, r y = r xy sin <fr with fixed 
r z = a\R\, |i ?| 2 = r 2 y +r\. Averaging over a uniform 
distribution for </> gives an unnormalized propagator, 


p' = p + a sin(2| R\) 


c y T ic x 


Cy ICx 


( 21 ) 


+ sin 2 |J?| 


2 (a 2 - l)c 2 


-(a 2 + l)(ca; - icy) 


-(a 2 + 1 )(c x + icy) -2(a 2 - l)c 2 
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FIG. 6 . A schematic drawing of a molecular beam experiment where an entangled state (axspin) is subjected to a random 
perturbing field before being detected at angle 9. 


When a = 1, r xy = 0, so there are no random forces. In 
this case, the density matrix above reduces to the exact 
equation of motion for a particle precessing in a magnetic 
field along the z-direction. 

The temperature-dependence of the stationary state 
agrees with the qualitative results on the harmonic oscil¬ 
lator from Sec. lIV Bll It resembles e ~ l3H , while deviating 
toward larger effective temperatures for large (3 when the 
noise is strong (a « 1 , not shown). 

Extracting the matrix elements from Eq. [22] leads to 
the transition tensor, 


1 — a 

1 -b 

aeP E 


1 - b* 

ae~P E 

1 — a 


a = (1 — a 2 ) sin 2 |i?| 
b = (1 + a 2 ) sin 2 | R\ — ia sin(2|I?|) 

The outer indices are ij' for p' , with each 2x2 sub-block 
indexed by i,j for p. 

The expectation of the energy change per unit time 
can be found from ( dH ) = — d Io f ■ For a state with 
P(|) = P{ t) = random collisions cause the energy to 
decrease on average, 

{dH) = -^log(l - a + acosh(/3E)) 

—2 aE sinh(/3P) 

1 — a + acosh(/3E) ’ 

The process of collapse happens after even a single 
time-step because the modeler’s state of knowledge is now 
described by a weighted set of possible wavefunctions, 
whereas the actual ensuing dynamics will be associated 
with only one of those choices. 

Figure [7] shows the time-course of work and system en¬ 
tropy for the dynamics of a spin system starting from the 
x superposition state (left panel), p = [ \ j;] /2, and from 
a statistical mixture of up and down states (right panel), 
p = [ J i ] /2. Energy and information entropy differences, 
dH and dSi n f at each time step are differences between 
thermostatic state functions, Tr[ffp(t)] and Tr[—logp]. 
The total entropy increase ( dS tot , Eq.[l2j) is a measure of 
irreversibility, not a state function. Due to its definition, 
it contains an exact contribution, dSi n {. Subtracting this 


gives a quantity which we term dQ , to represent the heat 
given off to the environment. 

From the perspective of z as the unique axis, both 
starting states look identical - giving a 50% probabil¬ 
ity of detecting the starting state as ‘up.’ However, the 
information entropy of the first matrix is the minimum 
possible (zero), while the second is the maximum possible 
(log 2 ). As time progresses in the left panel, the entropy 
of the state first increases due to random collisions be¬ 
fore finally decreasing to settle down into an equilibrium 
state. However, the energy continually decreases as colli¬ 
sions tend to rotate the spin into the z direction on aver¬ 
age. This gives an initially negative ‘heat,’ since knowl¬ 
edge of the pure starting state could have been used to 
cool a properly constructed external environment, after 
the cross-over region, the density matrix is compressed 
toward the up-state, decreasing the information entropy 
and requiring a release of heat into the environment. 

The initial heating region is completely absent when 
starting from the uniform distribution in the right panel 
of Fig. [3 Since the system starts in the maximally uncer¬ 
tain state, it continually gives off heat while lowering its 
energy. Initially, the compression (measured by decrease 
in information entropy) is slow with large decreases in 
energy. Eventually, the system settles down into a final 
compression phase similar to the left panel. 

The dotted lines are plotted using the path relation for 
the total entropy derived in Eq. [14] It can be seen that 
the total entropy increase per unit time is always positive, 
reflecting the trade-off between work and entropy. Work 
has to be done to compress the distribution toward the 
up-state. In the absence of environmental noise, super¬ 
position makes it impossible for this to happen. Here, it 
occurs because our probability assignment for collisions 
accounts for the fact that the energy of the system plus its 
environment is conserved. Because of this, the environ¬ 
ment is not allowed to impart arbitrarily large amounts 
of energy to the system. The total entropy has to in¬ 
crease when work is done in this ‘thermal’ way, since the 
number of allowable collision events decreases with time. 


V. DISCUSSION 

It is at first surprising how closely the dynamics of 
the energy-weighted propagator matches the model with 
an explicit cavity and measurement process. Both show 
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FIG. 7. Entropy and heat production per unit time from the spin system (Fig. [6} during wavefunction collapse. The left panel 
starts from the pure aj-state, while the right panel starts from the uniform directional distribution (p = I 2 ). The ‘heat released’ 
is defined as the difference, dQ = dStot — dSUf. Parameters: The oscillator energy was ftw = 1, with ft = 2. The distribution 
for r in G° was integrated over the full 3D normal distribution using Gauss-Hermite quadrature and D = 10 -3 . 


a nearly exponential decay toward an eventual steady 
state. The steady-state is very nearly a Boltzmann dis¬ 
tribution at high temperatures or weak coupling. This 
result should be expected to hold independently from 
the specifics of the environmental noise process. The 
nonequilibrium dynamical method here thus fulfills its 
goal of providing a first-principles model for the dynamic 
energy exchange processes which derives the equilibrium 
distribution. Both models also show deviations toward 
higher effective temperatures due to environmental cou¬ 
pling - an effect observed in experiments, flij 

The low-temperature behavior of the stationary state 
is a physically real effect. It occurs here because of the 
continual random influence of the environment, remov¬ 
ing the possibility of approaching absolute zero temper¬ 
ature. It will have an imperceptible effect on the energy- 
temperature equation of state, since the behavior of the 
system is dominated by the ground state at these low 
temperatures. Unless the interaction is unreasonably 
large, this method still gives the probability of being in 
the ground state at 1 to ~7 decimal places. 


A. Relation to Entropy and Work Fluctuation 
Theorems 

Both the definition (Eq. [Tdl) and its interpretation dif¬ 
fer from the usual fluctuation theorems. [Tl], [33 Conven¬ 
tional fluctuation theorems are defined with respect to a 
Boltzmann-like stationary distribution and an exact, re¬ 
versible, dynamical process (or a stochastic process obey¬ 
ing detailed balance for that stationary distribution). 
The transition entropy approach pursued here begins by 
postulating a stochastic transition process and only then 


determines its equilibrium and nonequilibriumn proper¬ 
ties. 

Because of this difference in origin, the entropy defi¬ 
nition also differs from other expressions related to the 
second law. m m m m The origin of the usual fluc¬ 
tuation theorem in exact dynamics leads to vanishing 
irreversible heat dissipation in the log-average, 

dStot = (logP(AE)/P(-AE)) 

= fd (AE — AA) = AS inf , 

where A is the Helmholtz free energy of a thermostatic 
equilibrium state at temperature ks/P allowed by the 
dynamics. Of course, the equilibrium relation between 
A and (E) = U sets dStot = ASi n f, and there is no irre¬ 
versible entropy production. In order to make statements 
about irreversibility, exact fluctuation relations have to 
be supplemented by a model for the thermalization of 
the environment. It is the entropy increase of the envi¬ 
ronment that is ultimately responsible for irreversibility 
derived in this way. 

Several work fluctuation theorems have been adapted 
to the quantum mechanical setting using moment gener¬ 
ating functions, 

Tr [e~ ilH U{t)e ilH pU\t)} , 

for the energy difference between initial and final 
states. Hi These are closely related to the energy- 
weighted propagator derived here (Eq.[12]). They are the 
natural analogues to the definition of work distributions 
used for deriving the classical fluctuation theorems. The 
twist of perspective in this work shows that their success 
can also be viewed as justification for postulating a statis¬ 
tical law of random transitions. Viewed in this way, the 
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exponential averages are not just a calculational device 
for comparing to some ‘instantaneous reference distribu¬ 
tion’ but are also useful statistical likelihoods for observ¬ 
ing transitions in energy space. The results of this work 
show that these likelihoods can be relied upon to generate 
a physical steady-state distribution for an open quantum 
system - without any reference to the Boltzmann-Gibbs 
equilibrium. 

VI. CONCLUSIONS 

This work provides a new perspective on the dynam¬ 
ics of open quantum systems. Three complete examples 
were analyzed. Understanding this perspective motivates 
a new statement of the central question for nonequilib- 
rium systems: “What is the most logically consistent way 
to represent an experimenter’s subjective probability of 
the system-environment interaction when only the sys¬ 
tem’s trajectory is known?” Utilizing a model of the 
maximum entropy type for this interaction gives an im¬ 
plicit representation of the environment based only on its 
role as a source of noise and dissipation. In addition, the 
dynamics is qualitatively very similar to results with a 
more explicit representation of the environment. 

A second question was also considered, related to the 
most complete description of the system in-between in¬ 
teractions with the environment. The firmly established 
measurement hypothesis of quantum theory was used to 
answer it. A single wavefunction is a complete descrip¬ 
tion of an isolated system, but a density matrix repre¬ 
sents the maximal combination of physical and statisti¬ 
cal information available for measurement. Agreement 
with the peculiar property that individual wavefunctions 
do not represent mutually exclusive events was used to 
select an admissible form for the propagator. The con¬ 
clusion was to constrain the energy change induced by 
the operator rather than the energy change of the wave- 
function itself. Subtleties with normalization of nonlinear 
models based on the density matrix lead to the final form 
for Eq. [HI 

A completely linear propagator, G y N (Eq. [TTJ) was 
mentioned as a possibility, but rejected because it did 
not clearly connect initial and final wavefunction states 
between time-steps. Although its stationary state is sim¬ 
pler to find numerically, preliminary results (not shown) 
display a minimum effective temperature. At very low 
temperature (high A), the stationary state probabilities 
start to distribute uniformly across energy levels. 

The dynamics of the favored propagator, Gn lEq.fTSl). 
were compared to a reduced model for thermalization of 
a trapped atom in an optical cavity. In addition to show¬ 
ing similar relaxation rates and dynamical behavior, the 
energy-weighted propagator showed mathematical prop¬ 
erties consistent with maximum entropy thermodynam¬ 
ics. Specifically, the transition free energy (Eq. [13j acts 


as a generating functional for the distribution of energy 
flux (Fig. 0. 

In addition, the transition entropy (Eq. [nj describes 
the trade-off between cooling of the system (lowering its 
information entropy) and heat release into the environ¬ 
ment. Running a collapse experiment on a spin in a pure, 
superposition state initially absorbs heat from the envi¬ 
ronment, whereas the same experiment on a completely 
random spin continually decreases entropy, requiring a 
continuous heat release into the environment (Fig. 0. 
The physical manifestation of this difference is because 
incoming random light will be scattered in a preferen¬ 
tial set of directions by the known spin. The transition 
entropy serves as a clue to this opportunity. However, 
actually using this initial absorption of heat requires a 
way to make use of the fact that the particle begins in a 
known state, e.g. by capturing the scattered light. Oth¬ 
erwise, the known force exerted by the particle is silently 
converted into heat. Again, we find that entropy is a 
consequence of the interplay between the known and un¬ 
known. 

This work presents a consistent picture of the role of 
the observer and the environment throughout the process 
of wavefunction collapse. Collapse is not possible with¬ 
out a noisy environment. Interacting with that environ¬ 
ment, however, will cause further entanglement rather 
than collapse. Nevertheless, the interaction will remove 
coherences present in the system, while favoring ending 
states that are lower in energy. At this point, the ex¬ 
perimenter’s state of knowledge about the system is de¬ 
scribed by a density matrix that accounts for these dif¬ 
ferences in energy. However, the actual dynamics of the 
system will be due only to a single wavefunction belong¬ 
ing to that density matrix. The most logically consistent 
way to guess the system’s state after some amount of 
interaction has taken place is to choose a wavefunction 
from the density matrix using the usual quantum mea¬ 
surement rule. Put another way, wavefunction collapse 
never occurs, but the experimenter’s state of knowledge 
becomes heavily weighted toward individual wavefunc¬ 
tions representing a collapsed state. In this interpreta¬ 
tion, Schrodinger cat states are incomplete until some 
environmental interaction takes place. 

The reduction of the environment to a statistical noise 
source with a bare minimum number of properties is crit¬ 
ical to achieving this picture. Rather than starting from 
an exact dynamical process for the environment, it is 
possible to present a set of interaction Hamiltonians and 
derive a weighted propagator. It is immediately clear 
how to add time- and spatially-dependent models for the 
environment (e.g. correlated noise sources). Other max¬ 
imum entropy-type constraints, e.g. the change in the 
total momentum, y] . p,;, can also be included in the bias¬ 
ing exponential along with the energy weight. The result 
is a greatly simplified model of the essential elements for 
nonequilibrium quantum systems. 
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APPENDIX 

VII. COMPARISON WITH THE DENSITY 
MATRIX MASTER EQUATION 

This section shows that another well-known formal¬ 
ism for simulating the dynamics of open quantum sys¬ 
tems can be derived as a low-order approximation of the 
energy-weighted propagator (Eq. [12]). 

Master equations for the density matrix are the most 


popular means for describing the role of environmental 
noise in localizing particles on a microscopic scale. [50| 
Caldeira and Leggett began with a description of col¬ 
lapse of a density toward the diagonal induced by cou¬ 
pling to a bath of harmonic oscillators. However, this 
coupling alone does not provide a suitable mechanism 
for energy dissipation. Later, these decoherence models 
were extended to add dissipation by imposing an energy 
distribution on the bath variables. [5lL |52| 

The particular high-temperature limit derived by 
Caldeira and Leggett for external forces linear in x 
is. [lC 


! rrr , *7rr- i i 2 my . 1 , 2 .. 

P= + -^[xpx- -{x ,p\), 

(23) 


where the damping rate should be of dimension 7 ~ 
\J c /2 m. 

It will be shown that the terms of Eq. [23] are in 
1:1 correspondence with those appearing in the high- 
temperature limit of G y N [p\. In this limit, Gn and 
are identical. 

The first step is approximating the energy weighting 
on U r . Since the random forces, r, are of order yjAt, 
terms up to r 2 contribute, 


e -Aff/2 e -ifi 
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When applying G\p], the sign of A is switched on the right 
of p. This is denoted below by W, so 

GP\p] = J P°(r)W(U r )Z~ 1 / 2 pZ~ 1 / 2 W(U r ) dr. (26) 


The r th contribution to G\p] becomes, for small At and 
dropping terms linear in r (which vanishes on averaging), 
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(27) 

(28) 


Of course, W(I) = I and W(H) = H. 

Taking the high-temperature limit allows a first-order 
expansion of the damped position operator, 

W(x) ~ x - ^\H,x\ = x+ ^-p (29) 

Similarly, W(x 2 ) ~ x 2 + %£{x,p}. 
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The expansion of Z = f W(U/)W(U r )Po(r)dr is just 
as for G\p] above. Defining the variance as f r 2 Po(r)dr = 
a 2 At, 


Completing the calculation above gives 
2 

p=-^[H,p\ + ^(xp x -±{x 2 ,p} (31) 

+ ^(PP X ~ X PP ~ P})) + ^ p - 

Rearranging \{x,p},p] = 2(xpp — ppx — ihp) cancels the 
constant. The expression multiplying iX, responsible for 
dissipation, becomes [{p, p},x], yielding a final expression 
identical to Eq. [23] The identifications can be made, 


cr 2 = 2mq/P (32) 

7 = Xa 2 /2m (33) 


Z ~ / + ( L£tw{x)W(x) - ^¥(W(x 2 ) + W{x 2 )) 


= 1 
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This model therefore corresponds to a high-temperature 
limit of Eq. [9l where A = /?. Numerical results in this 
work confirm that the steady state for Eq. [9] at low tem¬ 
perature is A = P/2. The temperature discrepancy ap¬ 
pears to be due to dropping higher-order terms in /?. 






